Out-of-equilibrium phase re-entrance(s) in long-range interacting systems 
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Systems with long-range interactions display a short-time relaxation towards Quasi Stationary 
States (QSSs) whose lifetime increases with system size. The application of Lynden-Bell's theory 
of "violent relaxation" to the Hamiltonian Mean Field model leads to the prediction of out-of- 
equilibrium first and second order phase transitions between homogeneous (zero magnetization) 
and inhomogeneous (non-zero magnetization) QSSs, as well as an interesting phenomenon of phase 
re-entrances. We compare these theoretical predictions with direct A''-body numerical simulations. 
We confirm the existence of phase re-entrance in the typical parameter range predicted from Lynden- 
Bell's theory, but also show that the picture is more complicated than initially thought. In particular, 
we exhibit the existence of secondary re-entrant phases: we find un-magnetized states in the theo- 
retically magnetized region as well as persisting magnetized states in the theoretically unmagnetized 
region. 

PACS numbers: 



I. INTRODUCTION 

In statistical physics, phase re-entrance is a quite typ- 
ical phenomenon occurring in many physical systems, 
such as spin-glasses, colloids and polymers, in which 
there is a competition between different entropic terms 
[U, 0, H, 0, H, Hi- A phase re-entrance is normally as- 
sociated with inverse melting, a counterintuitive phe- 
nomenon in which isobaric addition of heat causes a dis- 
ordered (e.g., liquid) phase to crystallize, the reverse of 
the usual situation. Phase re-entrance occurs when, pro- 
viding additional heating to the system, the latter under- 
goes a new transition, from the ordered to the disordered 
phase. The phenomenon of phase re-entrance has been 
widely studied at thermodynamic equilibrium in systems 
whose constituents interact through short-range forces. 

In this paper we give evidence to the existence of phase 
re-entrance also in the case of long-range interacting sys- 
tems in out-of-cquilibrium dynamical conditions. 

Long-range interactions are such that the two-body 
interaction potential decays at large distances with a 
power-law exponent which is smaller than the space di- 
mension. The dynamical and thermodynamical proper- 
ties of these systems were poorly understood until a few 
years ago, and their study was essentially restricted to 
astrophysics (stellar systems) and two-dimensional tur- 
bulence (large-scale vortices) Q- Later, it was recog- 



nized that long-range systems exhibit universal, albeit 
unconventional, equilibrium and out-of-equilibrium fea- 
tures Q. It is for instance well known that such sys- 
tems get t rapped in long-lasting Quasi-Stationary States 
Ea, ini, [H, El, Q [li, [ill, before relaxing to 



(QSS) 

thermal equilibrium. The duration of a QSS increases 
with the number of particles N in the system. Remark- 
ably, when the thermodynamic limit {N oo) is per- 
formed before the infinite time limit (t —f -l-oo), the sys- 
tem remains permanently trapped in QSSs. As a con- 
sequence, QSSs represent the only accessible experimen- 
tal dynamical regimes for systems composed by a large 
number of long-range interacting particles. This includes 
systems of paramount importance, such as non-neutral 
plasmas confined by a strong magnetic field Il7l.ll8ll , free- 
electron lasers [l^ and ion particle beams [20[. The 
ubiquity of QSSs has originated an intense debate [2l| 
about the mechanisms responsible for their emergence, 
their persistence, and their eventual evolution towards 
statistical equilibrium. In fact, QSSs keep memory of 
the initial condition and, as a consequence, they cannot 
be interpreted by making use of the classical Boltzmann- 
Gibbs approach. 

In a series of recent papers [22|, [H, [H, [H, , an ap- 
proximate analytical theory based on the Vlasoy equation 
and inspired by a seminal work of Lynden-Bell [27| in as- 
trophysics has been proposed. This is a fully predictive 
approach, enabling one to explain the emergence and the 
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FIG. 1: Magnetization as a function of time, for systems 
with different sizes. The continuous, dashed, short-dashed, 
dotted and dot-dashed lines correspond, respectively, to TV = 
2 • 10^, 10*, 5 • 10^ 2 • 10^ 10^ We see that QSSs start approx- 
imately at t = 40, and their duration increases with the the 
system size; eventually, they relax to Boltzmann equilibrium 
(thick line). Simulations are performed starting from a two- 
level distribution with energy U — 0.6400, and averaging over 
different system realizations with the same initial distribution 
(50, 100, 200, 500 and 1000, respectively). 



properties of QSSs from first principles [28| . 

In this paper we utilize a well-known hamiltonian toy 
model, the so-called Hamiltonian Mean Field (HMF) 
model pgj . to demonstrate that phase re-entrance(s) may 
also occur in a long-range interacting system, dynami- 
cally trapped in a QSS. The HMF describes the motion 
of N rotators, coupled through an equal strength cosine 
interaction. The Hamiltonian reads 

N N 

^=oE^'' + ^E[i-^°«(^.-^o] (1) 



where 6j represents the orientation of the j-th rotator 
and Pj stands for its conjugated momentum. To monitor 
the evolution of the systems it is customary to introduce 
the magnetization, an order parameter defined as 



M 



N 



where nii = (cos 6*^, sin 0.;). (2) 



The infinite-range coupling between rotators, provides 
the system with all typical characteristics of a long-range 
system, as clearly displayed in Fig. [T] Here, the magneti- 
zation is monitored as a function of time: after an initial 
"violent" relaxation, the system reaches a QSS, which 
is followed by a slow relaxation towards Boltzmann sta- 
tistical equilibrium. The larger the system, the longer 
the intermediate phase where it remains confined before 
reaching the final equilibrium. 

The paper is organized as follows. In Sec. [TTl 
we present the continuous Vlasov picture and review 
the maximum entropy principle based on the Lynden- 



Bell approach. This theoretical setting is used to ob- 
tain the HMF phase diagrams in different representa- 
tions from which out-of-equilibrium phase transitions 
m, m m, m, ill and phase re-entrance [H, ^ can 
be predicted. In Sec. IIIIl these theoretical predictions 
are compared to A^-body simulations based on ([T]). Fi- 
nally, in Sec. IIVI we sum up our results and draw our 
conclusions. 



II. OUT-OF-EQUILIBRIUM PHASE 
RE-ENTRANCE: THE PREDICTION OF 
LYNDEN-BELL THEORY 

A. General theory and two-levels approximation 



In a recent series of papers [22|, |23|, l2J, |25|, l26| , an ap- 
proximate analytical theory based on the Vlasov equa- 
tion has been proposed for the HMF model stemming 
from the seminal work of Lynden-Bell [l^]. This is a 
fully predictive approach, justified from first principles, 
which captures most of the peculiar traits of the HMF 
out-of-equilibrium dynamics. The philosophy of the pro- 
posed approach is reviewed in the following. 

In the limit of oo, the HMF dynamics can be 

formally described using the Vlasov equation 



%+P%~ iM4f] sine - My[f] cose) ^ 



0, (3) 



where f{9,p^t) is the one-body microscopic distribution 
function (DF), and the two components of the magneti- 
zation are respectively given by 



MM] 



f cos 9d9dp, 



(4) 



My[f] = J f sin ed9dp. 



The mean field energy can be expressed as 



U^^ J fp'dOdp- 



Ml + Ml 



(5) 



Working in this setting, it can be then hypothesized that 
QSSs correspond to stable stationary equilibria of the 
Vlasov equation on a coarse-grained scale. Lynden-BeU's 
idea goes as follows. The Vlasov dynamics induces a 
progressive filamentation of the initial single particle dis- 
tribution profile, i.e. the continuous counterpart of the 
discrete A^-body distribution, which proceeds at smaller 
and smaller scales without reaching an equilibrium. Con- 
versely, at a coarse-grained level, the process comes to an 
end, and the distribution function ./qssI^jPi averaged 
over a finite grid, eventually converges to an asymptotic 
form. Following Lynden-Bell, one can then associate a 
mixing entropy to this process. Assuming ergodicity (i.e., 
efficient mixing), /q55(0,p), is obtained by maximizing 
the mixing entropy, while imposing the conservation of 
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Vlasov dynamical invariants. It is worth emphasizing 
that the prediction of the QSS depends on the details 
of the initial condition [sO], not only on the values of 
energy and mass as for the Boltzmann statistical equi- 
librium state. This is due to the fact that the Vlasov 
equation admits an infinite number of invariants, i.e. the 
Casimirs or, equivalently, the moments A4n = J f^dOdp 
of the fine-grained distribution function. The first mo- 
ment Ml = M. is just the normalization of the distribu- 
tion function, and we can refer to it as the conservation 
of the total mass. 

In the case of a two-levels initial condition, where the 
fine-grained DF takes only two values / = /o and / = 0, 
the invariants reduce to M and [/, since the moments 
Ain>i can all be expressed in terms of M. and /o as 

Mn = jT'dedp = Jf^-^xfd0dp = Jf^^JdOdp = 

/q^^M. For the specific case at hand, the Lynden-Bell 
entropy is then explicitly constructed from the coarse- 
grained distribution function / and reads [l^, [3l| : 



s[f] = - dpde 



/" / 

-^In- 

JO JO 



fo 



fo 



We thus have to solve the maximization problem 



max{5[/] 
/ 



U[f] = U,M[f] = l}. 



(6) 



(7) 



This maximization problem assures that the distribu- 
tion function is thermodynamically stable (most prob- 
able macrostatc) in the sense of Lynden-Bell The 
maximization problem ([7]), for an arbitrary functional of 
the form S[f] ~ ~ J C{f) dOdp where C is convex, also 
forms a criterion of formal nonlinear dynamical stabil- 
ity with respect to the Vlasov equation |32|. In the two 
levels approximation, where the Lynden-Bell entropy is 
a functional of /, the criteria of dynamical and Lynden- 
Bcll thermodynamical stability coincide. From Eq. ([6]), 
we write the first order variations as 



5S ~ I35U - a6M = 0, 



(8) 



where the inverse temperature (3 = \/T and the "chem- 
ical potential" a are Lagrange multipliers associated to 
the conservation of energy and mass. Requiring that this 
functional is stationary, one obtains the following distri- 
bution (ii,!!!: 



fQSsiO,p) 



fo 



1 + ef^^P - ^-[/QSs] cose- My[f Qss] sin0) + a 

(9) 

As a general remark, it should be emphasized the 
"fermionic" form of the distribution, which arises be- 
cause of the form of the entropy. Notice also that the 
magnetization is related self-consistently to the distribu- 
tion function by Eq. ^ , and the problem hence amounts 
to solving an integro- differential equation. In doing so, 



we have also to make sure that the critical point corre- 
sponds to an entropy maximum, not to a minimum or a 
saddle point. Let us now insert expression ^ into the 
energy and normalization constraints and use the defini- 
tion of magnetization (|4]). F urther, defining A = e" and 
m= (cos 6*, sin 0) yields [3^ : 



fo 
fo 
fo 
fo 



/ d6l/_ 



3/2 



1/2 

2 \P 

(3 



1/2 



(Ae 



-,3M-m\ 



(10) 



d0/i/2 (Ae-^^-) - 



U 



AP - 1 



- / d6'cos6i/_ 



- / d6lsin6'/_ 



1/2 



1/2 



(Ae 
(Ae 



-/3Mm 



-/3M-m 



) = M„ 

) = My, 



where we have defined the Fermi integrals 



In{t) 



X" 



1 + ie^ 



-dx. 



(11) 



Its asymptotic limits are recalled in |23 | . 

If we consider spatially homogeneous configurations 
{Mqss ~ 0); the Lynden-Bell distribution becomes 



fo 



/QSS(P) = 



/3p2 /2 ■ 



(12) 



In the non degenerate limit A —>■ +oo, the latter reduces 
to the Boltzmann distribution / = (P/2ny^^e~^P 1"^ and 
in the completely degenerate limit A ^ 0, it becomes a 
step function: / = /o for p < l/(47r/o) and / = other- 
wise. If wc make use of Eqs. (fTO)) . we get the following 
parametric relation, with parameter A, between the in- 
verse temperature /3 and the energy I] (for a fixed value 
of /o) [21: 



U - 



A/2 (A) 

^-l/2(A)' 

f3 



G(A), 



(13) 



1/2 



(A)^ 



where G{X) is a universal function monotonically increas- 
ing with A. A solution of the above equation certainly 
exists provided that 



(14) 



To compute G{0) one can use the asymptotic expansions 
of the Fermi integrals. This yields G{0) = 1/12. There- 
fore, the homogeneous Lynden-Bell distribution with 
fixed fo exists only for [2^ : 



U > f/™„(/o) 



1 



96^2/2 



The completely degenerate limit (A 
geneous distribution corresponds to: 



(15) 

0) of the homo- 
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(i) U = Uminifo) and (3 ^ +00 for any /q. 
The non degenerate limit (A +00) of the homoge- 
neous distribution corresponds to: 

(i) U — *■ +00 and /3 — > for any /q. In this hmit, the 
caloric curve ([13]) takes the classical form U ~ 1/(2/3) + 

1/2 [13, Hi]. 

(ii) fo — > +00 for any U. 

Let us now address the problem of stability of the ho- 
mogeneous Lynden-Bell distribution. In , it has been 
shown that the homogeneous Lynden-Bell distribution is 
stable if, and only if: 

/_i/2(A)A|r i/2(A)| < j^^. (16) 

If the distribution function satisfies p6|) then it is both 
Lynden-Bell thcrmodynamically stable (entropy maxi- 
mum) and formally non-linearly dynamically stable [33 | . 
Otherwise, it is Lynden-Bell thcrmodynamically unsta- 
ble (saddle point of entropy) and linearly dynamically 
unstable. For a given fo, the relation (|16p determines 
the range of A for which the homogeneous distribution 
is stable/unstable. Then, using Eqs. p3)) . we can deter- 
mine the range of corresponding energies. Specifically, 
the critical curve Uc{fa) separating stable and unstable 
homogeneous Lynden-Bell distributions is given by the 
parametric equations [2^ : 

/_i/2(A)A|r,/2(A)| = (17) 

1 _ 1 -^1/2 (A) 

2 8^Vo'^-i/2(A)3' 

where A goes from (completely degenerate) to -l-oo 
(non degenerate). In fact, the criterion only proves 
that / is a local entropy maximum at fixed mass and 
energy. If several local entropy maxima are found (for 
example, homogeneous and inhomogeneous Lynden-Bell 
distributions), we must compare their entropies to de- 
termine the stable state (global entropy maximum) and 
the metastable states (secondary entropy maxima). For 
systems with long-range interactions, metastable states 
have in general very long lifetimes, scaling like e^, so 
that they are stable in practice and must absolutely 
be taken into account [sj, [s^l- For this reason, (out- 
of-equilibrium) stability diagrams do not coincide with 
phase diagrams. In fact, the latter require a careful in- 
vestigation of metastable states. 

B. Phase diagram in the {fo,U) plane 

The phase diagram of the Lynden-Bell distribution in 
the {fo,U) plane is shown in Fig. [21 We have also 
plotted the stabihty curve Udfo) of the homogeneous 
phase (split in two parts, ?7c(/o) and U"{fo)) defined by 
Eqs. p7)) and parameterized by A. On the left of this 
curve, the homogeneous phase is stable (maximum en- 
tropy state) and on the right of this curve it is unstable 



(saddle point of entropy) [s^. For /o — > -l-oo, we are in 
the non-degenerate limit A — > -l-oo and the stability cri- 
terion HH) is equivalent to A^ > A-K^f^. Using Eq. p^ . 
this yields U > Uc = 3/4. This is the critical energy 
associated to the Maxwell distribution. On the line of 
minimum energy U ~ Uminifo), we are in the completely 
degenerate limit (A 0) and the stability criterion ([TBI) 
is equivalent to fo < (/o)c = l/(27rV2). Using Eq. ([H)) . 
this yields U > Uc = 7/12. This is the critical energy 
associated to the spatially homogeneous water-bag dis- 
tribution (step function). Therefore, the minimum en- 
ergy curve Uminifo) crosses the stability curve Udfo) at 
((/o)c,i7c)^ (0.1125,0.5833). 

If we now take into account Lyndcn-BeU's inhomo- 
geneous states, solving numerically Eqs. (jlOp . we find 
that the phase diagram displays first-order and second- 
order phase transitions. The curve Udfo) splits in two 
curves U'dfo) and J7"(/o). In the case of a second-order 
phase transition, the stability threshold corresponds to 
the transition between a homogeneous and an inhomo- 
geneous distribution. The second-order phase transition 
corresponds to the branch ?7"(/o)- On the other hand, 
for a first-order phase transition, as we have the coex- 
istence of two entropy maxima, the stability condition 
of the homogenous phase is no more sufficient to find 
the transition line, which has to be calculated by making 
a comparison between the two entropy maxima. This 
procedure has been followed to plot the line Udfo) in 
Fig. [21 This line of first-order phase transition is reached 
when the homogeneous and inhomogeneous phases have 
the same entropy. The line Urifo) (first order) and the 
line U'^ifo) (second order) merge together at a tricriti- 
cal point, located at ((/o)*,?7*) ^ (0.10947,0.608) and 
corresponding to A* = 0.024. We have also plotted the 
curves Udfo) and Umetaifo) giving the lateral edges of 
the metastability regions for the homogeneous and inho- 
mogeneous phases (see figure caption for more details). 

In conclusion, the second order phase transition oc- 
curs for a range of values of C^(./o) bounded by the tri- 
critical point (C/*,(/o)*), and by Uc = 3/4, reached for 
/o +00. For U > Uc = 3/4, the Lynden-Bell theory al- 
ways predicts a homogeneous phase (for any value of fo)- 
For fo < ifo)*, the homogeneous phase is always stable 
(for any U > Uminifo))- For fo > (/o)c, the homoge- 
neous phase is unstable for Uminifo) ^U < Udfo) and 
stable for U > Udfo)- The first order phase transition 
occurs for a range of Udfo) bounded by the tricritical 
point (C/*,(/o)*) and by the point iifo)r,Udifo)r)) - 
(0.1098,0.5875). As can be seen in Fig. [2 the theory 
predicts a phase re-entrance, for a set of values of fo 
G [(/o)*, ifo)c\- This means that, increasing U in the di- 
agram at fixed fo G [(/o)*, (/o)c], the homogeneous phase 
is (meta)stable for Uminifo) < U < Udfo), unstable for 
Udfo) <U < U'^ifo), and stable again for U > C/"(/o)- 
Note that in the metastability region Uminifo) < U < 
Udfo) J the system can be found either in the homoge- 
neous or inhomogeneous phase depending on how it has 
been prepared initially (recall that metastable states are 
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FIG. 2: Phase diagram in the (/o, U) plane. The homoge- 
neous phase only exists above the line ?7mm(/o). The sta- 
bility curve Ucifo) is parameterized by A. For A — » (com- 
pletely degenerate limit), we get fo — (/o)c = l/(27r\/2) and 
Uc = 7/12. For A — > +oo (non degenerate limit), we get 
fo —> +0O and Uc = 3/4. On the left of this curve, the ho- 
mogeneous phase is stable and on the right of this curve it 
is unstable. The stability curve is divided in two parts, i.e. 
Uc{fo) and [/"(/o), by the tricritical point (full round dot) lo- 
cated at ((/o)t,i7*). The continuous line corresponds to the 
second-order transition line while the dotted lines correspond 
to the borders of the metastable region. The thick line rep- 
resents the first-order transition line. All these lines divide 
the diagram in four regions. In region (I), the homogeneous 
phase is stable and the inhomogeneous phase does not exist; 
in (II), the homogeneous phase is stable and the inhomoge- 
neous phase metastable; in (III), the homogeneous phase is 
metastable and the inhomogeneous phase stable; in (IV) the 
homogeneous phase is unstable and the inhomogeneous phase 
stable. Uminifo) is the line below which the homogeneous 
phase does not exist, and UMiN^fo) is the lowest accessible 
value of energy for a rectangular water bag initial condition 
(see Sec. IIIEp . The square dot is {{fo)c,Uc{{fo)c)), the dia- 
mond is {{fo)m.,Um) (Sec. IIIEp . and the empty round dot is 
((/o)r, f/r((/o)r)). For fo G [(./o)., {fo)c] there is a re-entrant 
phase. 



highly robust for systems with long-range interactions). 
By contrast, for U'^{fo) < U < U^'{ff)), the theory pre- 
dicts an inhomogeneous phase and for U > ?7"(/o) a 
homogeneous phase. 



Now, many numerical simulations of the TV-body sys- 
tem [23I , or of the Vlasov equation [23| , have been per- 
formed starting from a family of rectangular water-bag 
distributions. The latter correspond to assuming a con- 
stant value /o inside the phase-space domain D: 

D = {{e,p) e [-7r,7r] X [-00, +00] I \e\ < Ae, \p\ < Ap} 

(18) 

and zero outside. Here < Ad < tt and Ap > 0. The 
normalization condition results in 



./o = 



1 



iAOAp' 



(19) 



Notice that, for this specific choice, the initial magneti- 
zation Mq and the energy density U can be expressed as 
functions of A9 and Ap as 



U ■ 



[Apf , l-(Afo)^ 



sin(A6i) 
AO 



(20) 



(21) 



For the case under scrutiny, < A/q < 1 and U > 
Umin{Mq) = (1 - Af2)/2. The energy Umin{Mq) rep- 
resents the absolute minimum accessible energy for a 
rectangular water-bag distribution with magnetization 
Mq. The initial configuration is completely specified by 
the variables (A0, Ap) or, equivalently, by the variables 
(A/q, U). On the other hand, for the determination of the 
Lynden-Bell equilibrium state, only the variables (/o, U) 
matter. Now, we note that different values of {Mq,U) 
can correspond to the same (/o, U) and, consequently, to 
the same Lynden-Bell equifibrium (see Sec. IIIPp . There- 
fore, the use of these variables leads to some redundances. 
Nevertheless, their advantage is that they are more di- 
rectly related to physically accessible parameters. In any 
case, it is of interest to compare the two phase diagrams 
in (/o, U) and (Afo, U) planes to see their similarities and 
differences. 

For the rectangular water-bag initial condition, using 
Eqs. (jl9|) and ([20|) . we can express /o as a function of 
A/q and U by 



C. Phase diagram in the {Mo,U) plane 

The preceding results are valid for any initial condition 
with two phase levels f ~ and f ~ fo, whatever the 
number of patches and their shape. In the two-levels case, 
the relevant control parameters are {fo,U) They 
fully specify the Lynden-Bell equilibrium state from the 
initial condition. This means that, assuming ergodicity, 
the system remembers the initial condition through the 
values of these parameters. In this sense, the general 
phase diagram in the two-levels case is the one repre- 
sented in Fig. [21 



[(2[/- l)(A6i)2 



'Aoy 



(22) 



where A9 is related to Mq by Eq. (HI]). Inserting this 
expression in Eqs. (|13p . we obtain after some algebra the 
caloric curve T(U) for fixed Mo parameterized by A: 



U 



sin^ AO 



t2 I- 



? A0 



1/2W 
^2 



2(A6i)2 



2(A6') 



2 -^1/2 (A) 

^-l/2(A) 



(23) 
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Using Eqs. (fT5|) and p2l) . wc find that the homogeneous 
Lynden-Bell distribution with fixed Mq exists if and only 
if 



U > U,mn{Mo) = \ 



sin^ Ae* 
■2 - (A6i)2 



(24) 



Tliis result ean also be obtained from Eq. (|23p by con- 
sidering the limit A — > 0. 

The completely degenerate limit (A = 0) of the homo- 
geneous distribution corresponds to: 

(i) \J = Umin{Mo) and f3 — > +oo for any Mq. 

(ii) Mo = for any U. 

The non degenerate limit (A +oo) of the homoge- 
neous distribution corresponds to: 

(i) U +00 and /? ^ for any Mq. In that case 
[/ ~ 1/(2/3) + 1/2. 

(ii) Mo = 1 for any U. In that case /o +oo. 

On the other hand, regrouping all the preceding re- 
sults, the critical curve Uc{Mq) separating stable and 
unstable homogeneous Lynden-Bell distributions is given 
by the parametric equations 



i/2(A)A|/ii/2(A)| = 



1 



U 



1 I1/2W 
8^Vo'^-i/2(A)= 



[(2C/- l)(A6i)2 +sin^ AO] 



Mq 



sin(A6l) 
A9 '■ 



(25) 



(26) 



(27) 



(28) 



where A goes from (completely degenerate) to +00 (non 
degenerate) . 

The phase diagram of the Lynden-Bell distribution in 
the (Mo, U) plane is represented in Fig. [31 We have first 
plotted the minimum accessible energy of the homoge- 
neous phase ?7mm(Mo) defined by Eq. ([M)) . We have also 
plotted the stability curve J7c(Mo) of the homogeneous 
phase defined by Eqs. (PSjl - ipS]) and parameterized by A. 
Above this curve the homogeneous [Mqss = 0) phase 
is stable (maximum entropy state) and below this curve 
it is unstable (saddle point of entropy). For Mq = 1, 
we are in the non degenerate limit A +00 (because 
/o +00) and the critical energy is C/c = 3/4 (Maxwell 
distribution). For Mq = 0, we are in the completely de- 
generate limit A = and the critical energy is Uc = 7/12 
(spatially homogeneous water bag). 

If we now take into account Lynden-Bell's inhomo- 
geneous states, solving numerically Eqs. pU)) . we find 
that the phase diagram displays first-order and second- 
order phase transitions. The second order phase tran- 
sition corresponds to the branch U^'{Mq) and the first 
order phase transition to the branch J7r(Mo). These 
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FIG. 3: Phase diagram in the control parameter plane 
(Mo, U) for a rectangular water-bag initial profile. Umin[Mo) 
is the absolute minimum energy, and the homogeneous phase 
only exists above the line Umin{Mo). The stability curve 
Uc{Mo) is parameterized by A. For A — > (completely de- 
generate limit), we get Mo = and Uc = 7/12. For A — » -l-oo 
(non degenerate limit), we get Mo = 1 and Uc ~ 3/4. Above 
this curve, the homogeneous phase is stable and below this 
curve it is unstable. The full dot is the tricritical point. In 
the inset is showed the region of the first order phase transi- 
tion, indicated by the line Ur{Mo), connected to the second 
order phase transition line C/"(Mo) by the tricritical point. 
The dotted lines represent the borders of the metastability 
region. 



two lines merge together at the tricritical point, lo- 
cated at ((Mo)*,C/*) ~ (0.1757,0.608) corresponding to 
A9^ = 2.656... Using Eq. ([^ . we readily check that this 
tricritical point ((Mq)*, U^) corresponds to the tricritical 
point ((/o)*, U^) in the (/o, U) plane. We have also plot- 
ted in the inset the lateral edges U^{Mo) and Umeta{Mo) 
of the metastability region associated to the first order 
phase transition. 

Therefore, the phase diagrams in (/o, U) and (Mo, U) 
planes are fully consistent (see Sec. IIIDI for more de- 
tails) and both display first and second phase transitions. 
The correctness of the above analysis is assessed in [2^ 
where numerical simulations are performed for different 
values of the system size N. The transitions predicted in 
the realm of Lynden Bell's theory are indeed numerically 
observed, thus confirming the adequacy of the proposed 
interpretative scenario. Note, however, that the physics 
is different whether we vary the energy at fixed /o or at 
fixed Mq. In particular, there is a "re-entrant" phase 
when we vary the energy at fixed /o [2^ but there is no 
"re-entrant" phase when we vary the energy at fixed Mq 



D. Connection between the two phase diagrams 

To make the connection between the phase diagram 
{fo,U) of Sec. Ill Bl and the phase diagram {Mo,U) of 
Sec. Ill CI we can plot the iso-Mo lines in the {fQ,U) 
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phase diagram or the iso-/o lines in the {Mo,U) phase 
diagram. 

Let us first consider the iso-Afp hnes in the {fQ,U) 
phase diagram. If we fix the initial magnetization Mq, 
or equivalently if we fix the parameter A^, the relation 
between the energy U and /o is 



UAeifo) 



1 



sin Ae* 



6(4A6i/o)2 2 V A6I 
Therefore, the iso-A'/o lines are of the form 

A{Ae) 



(29) 



(30) 



• sin Ae ^ 
, A0 ) 



with A{Ae) = and B{M) - 

which are easily represented in the (fo,U) phase dia- 
gram (see Fig. H]). As an immediate consequence of 
this geometrical construction, we can recover the min- 
imum energy of the homogeneous phase for a fixed initial 
magnetization Mq (or Ad). Indeed, for a given A9, the 
homogeneous phase exists iff UAeifo) > Uminifo) leading 
to 

7r2 - (A0)2 



/o < (/o)Ae 



(31) 



487r2 sin^ AO 

This corresponds to U > UminiMo) — UAe{{fo)Ae) lead- 
ing to 



U > C/,„„(Mo) = 



sin^ A9 
7r2 - (A6l)2 



+ 1 



(32) 



which is identical to Eq. ([M)) . Figure 0] is in good agree- 
ment with the structure of the phase diagram in the 
(Mo, U) plane. Indeed, along an iso-Mo line, we find that 
for large energies U > Uc{Mo) the homogeneous phase is 
stable and for low energies U < Uc{Mo) the homoge- 
neous phase becomes unstable. In that case, there is no 
re-entrant phase. We also note that for < (A/q)*, 
the phase transition goes from second order to first or- 
der. This corresponds to the case where the iso-Afp line 
crosses the tricritical point. 

Remark: we see on Fig. |4] that different iso-A/o lines 
can cross each other. This means that different initial 
conditions {Mq,U) and {Mq,U) can correspond to the 
same {fo,U) hence to the same Lynden-Bell distribu- 
tion. In other words, in Lynden-Bell's theory the couples 
{Mq,U) and {Mq,U) are equivalent. There is therefore 
some redundance in using the variables {Mq , U) instead 
of the variables (/g, U). Note, however, that the Lynden- 
Bcll prediction does not always work so that, in case of 
incomplete relaxation, the couples (Afg,J7) and (A/q,[/) 
may lead to different QSS. However, this happens for 
/o > 0.11053 . . . , i.e. in a parameter range which is only 
marginally interesting for our analysis. 

Let us now consider the iso-/o lines in the {Mo,U) 
phase diagram. If we fix the phase level /o, the relation 
between the energy U and Mq, or equivalently A6, is 



Ufom 



1 



6(4A0/o)2 



sin Ae* 
A9 



(33) 
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FIG. 4: Iso-Afo lines in the {fo,U) phase diagram. This 
graphical construction allows one to make the connection be- 
tween the {fo,U) phase diagram of Fig. [2] and the {Mo,U) 
phase diagram of Fig. |3] We can vary the energy at fixed 
initial magnetization by following a dashed line. The inter- 
section between the dashed line and the curve Uminifo) de- 
termines the minimum energy Um.in{Mo) of the homogeneous 
phase. The intersection between the dashed line and the curve 
Uc{fo) determines the energy Uc{Mo) below which the homo- 
geneous phase becomes unstable. 



Therefore, the iso-/o lines are of the form 
Cifo) 



1 /sinA6l 

2 [ AO 



1 

2' 



(34) 



with C{fo) ~ e(4joy^ ' ^hich are easily represented in the 
{Mo,U) phase diagram (see Fig. [5|). Recall that Ad is 
related to Mq by Eq. (pi]) . FigureOis in good agreement 
with the structure of the phase diagram in the {fo,U) 
plane. In particular, we can see that for a set of values 
of /o € [(/o)*) (/o)c], it intersects the curve Uc{Mo) twice 
leading to re-entrant phases. We also note that the iso- 
/o lines cannot cross each other contrary to the iso-Afo 
lines. 



E. Determination of the absolute minimum energy 
in the {fo,U) plane for a rectangular water-bag 
initial condition 



We recall that homogeneous Lynden-Bell distributions 
exist only for U > Uminifa)- However, there can ex- 
ist inhomogeneous Lynden-Bell distributions for U < 
Uminifo)- Let us determine the minimum accessible en- 
ergy UminHo) when we start from a rectangular wa- 
terbag initial condition. For fixed /o the energy of the 
initial condition is a function of A9 (or initial magneti- 
zation Afo) given by 



Ufom 



1 



6(4A0/o)2 



sin Ae* 
A9 



(35) 
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FIG. 5: Stability diagram in the {Mo,U) plane, with iso-/o 
lines. The thick lines are the two parts of the stability curve, 
the dotted line is the iso-/o line with fo — 0.1096 and the 
dash-dotted line is the iso-/o line with fo = 0.1100. 



Wc thus have to determine the minimum of this function 
for 0< A9 <n. First of all, the condition C/}^^ {A9) = 
is equivalent to 



/o = 



1 



v/48 sin(A6l) (sin(A6i) - A6icos(A6')) 



(36) 



This function is represented in Fig. [6l For fo < 0.11053... 
there is no solution and for fo > 0.11053... there are 
two solutions A9i and A92 corresponding to one local 
minimum and one local maximum (see Fig. [7]). For fo = 
0.11053... , we have A9i = A92 = 2.2467... . Then, we 
find that the local minimum is the absolute minimum iff 
Uf„{A9i) < Ufyin). This is the case if fo > {fo)m = 
0.12135... corresponding to an energy Um = 0.57167... 
(see Fig. [7]). For fo < (/o)m, the absolute minimum 
correspond to A^^ = tt. 

In conclusion, for fo < {fo)nn we find that 
UminUo) = Ufgiji) so that 



UMINifo) — Uminifa) 



1 



967r2/o' 



(37) 



For fo > (/o)m, we find that UminUo) = Uf„{A9i) 
where A^i is the smallest root of Eq. ([55]) . Combin- 
ing these equations, we find that the absolute minimum 
energy UminHo) is given in parametric form by 



MIN 



1 



sin(2A6li) 

2Ae'i 



/o 



^48 sin(A6li) (sin(A6li) - A6ii cos(A6li)) ' 



(38) 



(39) 



with < A6'i < 1.85063... . For fo -> +oo (non de- 
generate limit), we get A^i and ?7a//jv(/o) ^ 0. 
This is the energy corresponding to an initial condition 
f{9,p,t = Q) = 6(j>)5{9). 
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FIG. 6: Graphical construction determining the solutions of 
the equation U'f^^{A9) = 0. 




FIG. 7: Energy of the initial condition as a function of A9 
for different values of fo. From top to bottom: fo — 0.10, 
/ = (/o)m = 0.12135, fo = 0.14. 



Remark: for fo > 0.11053, we confirm on Fig. [7] that 
there can exists several initial conditions with the same 
/o and U but a different initial magnetization Mq. They 
lead to the same Lynden-Bell prediction. 



III. NUMERICAL RESULTS 

To assess the correctness of the above theoretical pre- 
diction about the existence of a phase re-entrance, we 
have performed direct numerical simulations of the HMF 
model (dl) for finite A^. For that, we have chosen fo in 
the interval {fo)* < fo < (/o)c and ran simulations at 
different energies. Results, for fo = 0.1096, are shown 
in Fig. [51 where both the theoretical and numerical val- 
ues of magnetization at the QSS, Mqss, are plotted as 
a function of energy. 

Simulations (dashed line) confirm the existence of a 
regime of phase re-entrance. However, the agreement 
with theory (continuous line) is mainly qualitative, as 
there is a systematic shift between the two curves, al- 
though the magnetization value of the main bump is 
consistent with the one predicted from Lynden-Bell's ap- 
proach. Moreover, simulations show the existence of two 
zones of magnetization revival at both sides of the cen- 
tral magnetized region. If we move at fo = 0.1100 (Fig. 
[5]), we find that one of the two bumps has grown; this 
confirms that the structure of the phase diagram is more 
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FIG. 8: Magnetization value at QSS, Mqss, versus energy 
for fo = 0.1096. Comparison between theory (continuous line) 
and simulations (dashed line). Simulations, done with TV — 
10^, are performed using a symplectic integration algorithm, 
and averaging the magnetization over the time window 40 < 
t < 140, and over 50 different realizations. For this fo, Umin ~ 
0.5878, Ur = 0.5955, = 0.6026 and f/" = 0.6131. 
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FIG. 9: Magnetization value at QSS, Mqss, versus energy 
for fo = 0.1100. For this fo, U^in = 0.5872, = 0.5980 and 
Uc = 0.6177. We plotted it for different sizes of the system: 
the continuous line correspond to TV = 10^, the dashed to 
N = 5- and the dotted to TV = 10^ 



complex than predicted by the theory, as we find the 
existence of additional phase re-entrances. Simulations 
performed using different numbers of particles show that 
the magnetization values of the central magnetized re- 
gion and of the two bumps do not depend on the system 
size. Instead, as expected, the curve offset goes to zero 
when the system size is increased, see Fig. [HI 

A possible explanation for the discrepancies between 
theory and simulations can be found by considering that 
the energy range in which phase-rentrance is observed 
is quite narrow. In fact, the iso-/o lines, in the inter- 
val {fo)* < fo < {fo)c, are very close to the theoreti- 
cal phase-transition curve, see Fig. \S\ This means that 
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FIG. 10: Magnetization value, Mqss, versus energy for fo = 
0.1130 > (/o)c. Comparison between theory (continuous line) 
and simulations (dotted line). Umin{fa ~ 0.1130) = 0.5826 
and U"{fo = 0.1130) = 0.6325. 



any possible (i.e., even small) disagreement between the 
theoretical and the numerical one may easily lead: a) 
to a further numerical phase rc-entrancc, if the iso-/o 
line crosses the numerical phase-transitions curves with- 
out crossing the theoretical one; b) to a larger numerical 
value of Mqss, if the iso-/o separates from the numerical 
curve, while staying close to the theoretical one. 

We also compared theory and simulations for higher fo 
(> (/o)c)- As shown in Fig. [TOl here theoretical and nu- 
merical results are close. This is in agreement with what 
is reported in [2^. In Fig. [TT] we plotted our numeri- 
cal results for a fo lower than {fo)*- For fo = 0.1085, 
close to the critical line Uc{fo), we observe a magnetized 
phase although Lynden-Bell's approach predicts a non- 
magnetized phase. For lower values of fo, homogeneous 
QSS are observed in agreement with the theoretical pre- 
diction (data not shown). 

To provide a complete picture of the whole phase dia- 
gram, wc carried out simulations on a grid in the {fo,U) 
plane and plotted the numerically obtained values of 
Mqss in color scale, see Fig. [T^l At first order, we ob- 
serve a fair agreement between theory and simulations. 
In particular, the predicted re-entrant phase phenomenon 
is clearly observed. This can be considered as a success 
of the Lynden-Bell theory. We also note that the region 
(III) of the phase diagram appears to be non-magnetized. 
It corresponds therefore to a local Lynden-Bell entropy 
maximum. This confirms our claim about the robust- 
ness of metastablc states. Note, however, that starting 
from different initial conditions (with identical values of 
U and fo), we could have found that the QSSs in this 
region are magnetized. Indeed, in the metastability re- 
gion, the selection between between local (metastable) or 
global (state) entropy maxima depends on a complicated 
notion of basin of attraction. Furthermore, wc also find 
some unprcdicted phenomena, as the additional phase rc- 
cntrance, for U ~ 0.605 (sec also Fig. [5]) and a persisting 
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FIG. 11: Magnetization value, Mqss, versus energy for /o — 
0.1085 < (/o).. For /o < (/o)., the Lynden-Bell approach 
predicts that the QSS should be non-magnetized so there is 
a disagreement with theory when /o is close to the tricritical 
point. 
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FIG. 12: Stability diagram in the fo-U space with numeri- 
cally calculated mean magnetizations. The dashed line is the 
stability curve. The theoretical re-entrant phase is clearly vis- 
ible as well as the second (unexpected) re-entrant phase. In 
addition to this interesting new re-entrance phase, the other 
main discrepancy is the persisting magnetized phase for low 
fo- 



magnetized phase for low fo (also Fig. [TT|) . 

We also studied the order of phase transitions, by plot- 
ting the probability histogram of Mqss sampled with 
300 different realizations. We show the results for two 
of the transitions occurring in Fig. [5] In Fig. I13[ one 
can see that for the transition at U ~ 0.5980, distri- 
butions are characterized by a double peak, which is a 
clear signature of a first-order phase transition. For the 
one at [/ ~ 0.6230, see Fig. [T31 the distributions are 
instead characterized by a single peak, which validates 
the prediction of a second-order phase transition. The 
two transitions at the boundaries of the smaller phase 
re-entrance, occurring around U ^ 0.605, not predicted 
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FIG. 13: Probability distributions of Mqss for different U 
values at fo = 0.1100. Here in (a) U = 0.5970, in (b) U = 
0.5980, (c) U = 0.5990, in (d) U = 0.6000. 
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FIG. 14: Probability distributions of Mqss for different U 
values at fo = 0.1100. Here in (a) U = 0.6150, in (b) U = 
0.6190, (c) U = 0.6230, in (d) U = 0.6270. 



by the theory, are found to be of first (at low energy) and 
second (at high energy) order (data not shown). 

Finally, we compared the analytical and numerical 
caloric curves (3{U) for a given value of fo. In the sim- 
ulations, the temperature has been calculated from the 
usual expression 



1 



<p^ >= / d6dpfQss{0,p)p'' 



(40) 



We note that the "kinetic" temperature defined by Eq. 
(Uni) does not coincide with the Lagrange multiplier (3 as- 
sociated to the energy conservation in the Lynden-Bell 
distribution This is due to the fermionic nature of 
this distribution. Therefore, in order to make the com- 
parison between simulations and theory relevant, we have 
calculated the theoretical temperature from the mean 
square momentum ()40|) averaged with the Lyndcn-BcU 
distribution given by Eq. ([9]). The results are reported 
in Fig. [T31 In continuity with the results of Fig. [51 the 
range of energies where the inhomogeneous phase appears 
is shifted with respect to the theoretical prediction. As 
a further point, we also notice the presence of a region 
with negative specific heat, both in the numerical and an- 
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FIG. 15: Comparison between theoretical (continuous line) 
and numerical (dashed line) caloric curves, for fo = 0.1096. 

alytical curves. To the best of our knowledge, this is the 
first time negative specific heat is observed out of equi- 
hbrium. Surprisingly, this phenomenon is here observed 
in correspondence of a second order transition line. 

IV. CONCLUSION 

In this paper, we have confronted the predictions 
[13, SHIM Hi] of a theory based on L^den-BeU's sta- 
tistical mechanics of violent relaxation [23| to the results 
of numerical experiments. The application of Lyndcn- 
Bell's theory to the HMF model predicts a re-entrant 
phase in the (/o, U) plane [l^l and, indeed, we observe 
it. It occurs for a narrow range of parameters which 
would have been difficult to find without such a theo- 
retical prediction. In this sense, this is a great success 
of Lynden-Bell's approach. The theory also predicts the 
correct value of the magnetization in the inhomogeneous 
phase and the correct order of the phase transition. This 
is again remarkable because the phase diagram displays 
first and second order phase transitions in a very nar- 
row range of parameters {fo,U). All these predictions 
are confirmed by direct A^-body experiments. We have 
also numerically observed that metastable states (local 
Lynden-Bell entropy maxima) can be very robust, so that 
they are stable in practice. This is a specificity of systems 
with long-range interactions [13, [l^ • 



However, we have also found some discrepancies with 
respect to Lynden-Bell's theory. In particular, numeri- 
cal simulations have demonstrated the existence of sec- 
ond re-entrant phases: a band of un-magnetized states 
in the theoretically magnetized region, as well as persist- 
ing magnetized states in the theoretically un-magnetized 
region. As a matter of fact, there is a systematic shift 
of the transition line with respect to theory. We must 
emphasize, however, the very small selected region of pa- 
rameters in Figs. [5] and [T^ This gives the impression 
of a big discrepancy although the discrepancy is in fact 
very small. 

Therefore, from these numerical experiments, we can 
conclude that the Lynden-Bell statistical theory gives a 
fair first order description of QSSs in the HMF model. 
However, for some initial conditions, there can be more 
or less severe discrepancies with respect to the prediction. 
This is a well-known fact in stellar dynamics [l^ and 
vortex dynamics (isj to which this theory was initially 
applied (see a detailed discussion in [13]) ■ Discrepancies 
with the Lynden-Bell theory have also been reported for 
the HMF model in [23| and [13|- These discrepancies 
are usually the result of an incomplete relaxation [23|, 
i.e. a lack of efficient mixing in the system phase space. 
Indeed, the Lynden-Bell theory is based on a hypothesis 
of ergodicity and the prediction fails (by definition) if 
the evolution is not ergodic. A detailed understanding of 
incomplete violent relaxation is still lacking and appears 
to be very difficult 

Another cause of discrepancy may be related to 
the proximity of the numerically considered parameters 
[fo, U) to the critical line and to the tricritical point. In- 
deed, it is well-known in equilibrium statistical mechan- 
ics that strong fluctuations are present close to a critical 
point, and that the mean field results cease to be valid in 
the vicinity of a critical point [13] ■ In the present case, 
we are studying out-of-equilibrium phase transitions and 
it is not clear if we can directly extend equilibrium results 
to that situation. Nevertheless, it is not unreasonable to 
expect that the theoretical results may be altered close to 
the critical line and this is indeed what we observe numer- 
ically. Further away from the critical line (i.e. for larger 
or smaller values of /o), we find a very good agreement 
with the Lynden-Bell prediction (see also [131) ■ These 
different observations concerning the success or the fail- 
ure of the Lynden-Bell theory are consistent with the 
discussion given in [13|- 
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